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Abstract 

We study the nonequilibrium photon production in the quark-hadron phase 
transition, using the Friedberg-Lee type solitons as a working model for quark- 
hadron physics. We propose that to search for nonequilibrium photons in the 
direct photon measurements of heavy-ion collisions may be a characteristic 
test of the transition from the quark-gluon to hadronic phases. 
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Quark-confinement is an unresolved problem in QCD physics. Nevertheless, it is gen- 
erally believed that quarks are deconfined at high energies. One of the primary goals of 
relativistic heavy-ion-collision experiments such as BNL-AGS and CERN-SPS is to create 
a hot central region and a dense fragmentation region in which quark-gluon plasma (QGP) 
may be formed. One of the ways for testing the occurrence of QGP is to observe the spec- 
trum of the emitted photons. An advantage of measuring the direct photons is that photons 
do not suffer from strong final-state interactions as hadrons do, so they can be used to 
monitor the initial stages of the colliding heavy ions 0. 

In relativistic heavy-ion collisions, a large number of tt° and rj are produced by soft 
QCD processes. Their two-photon decays constitute the major source of photon emission. 
Other photon sources arise from radiative decays of other mesons and baryon resonances, as 
well as hadron scattering processes. Also, if QGP is formed, photons will be produced via 
quark- ant iquark annihilations or quark-gluon Compton scatterings. It is an experimental 
challenge to separate out the photon yield into a part arising from the decays of produced tt° 
and rj and the "single photon" part from other sources. For instance, the part from decays 
of mesons can be subtracted by reconstructing the distribution of the two-photon invariant 
mass of all photon pairs. If this can be done, then the excess single photon spectrum can be 
a probe of the properties of the nuclear matter or QGP, or even a discriminator between the 
two phases. However, recent studies have shown that the two single photon spectra are 
similar, and distinguishing between them might need further high-precision direct photon 
measurements or other complementary measurements such as dilepton production and J/ ip 
suppression. 

In this Letter, we give an attempt to utilize the nonequilibrium nature of the photon 
emission associated with the quark-hadron phase transition (QHPT) to test the formation of 
QGP. In essence, we calculate the particle production from the release of latent heat of the 
phase transition, adopting the formalism of nonequilibrium field theory that is well-suited 
to study the dynamics of nonequilibrium processes To model the QHPT, we use the 

Friedberg-Lee (FL) phenomenological nontopological solition model |§, which is simple and 
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adequate enough for our present considerations. In fact, FL model has been applied to fit 
the hadronic properties || , as well as to discuss the QHPT of the early Universe [ |T0| , |rTf . We 
believe that the particle production considered here is generic to all Landau-Ginzberg type 
models of QHPT. Our main result will be the photon spectrum produced from the QHPT, 
and the experimental signature will be briefly discussed. 

We begin with the simplest FL model whose Lagrangian is given by f§ 



where ipi represents the quark field, np = 2 is the number of light flavors, a is a real scalar 
field, / is a coupling constant, and the potential U(a) is given by 



where a, b, and c are positive parameters, and B is the bag constant. Note that U(a) has 
a local minimum at a = 0, and an absolute minimum at o = o c which denotes the true 
vacuum with U(a c ) = 0, separated by a potential barrier. A typical potential is shown in 



which are identified with hadrons. Essentially, inside a hadron is a perturbative vacuum 
with a ~ decorated with localized free quarks, whereas outside is the true vacuum. 

The simple picture of the production of nonequilibrium photons from the QHPT is the 
following. Suppose that a thermalized QGP is produced initially in a relativistic nucleus- 
nucleus collision, which has a perturbative vacuum with o = 0. It then expands and cools, 
and undergoes hadronization in a phase transition. Here we assume that the transition is 
weakly first-order as consistent with lattice QCD calculations As such, transition to 
a state with a ~ fa (see Fig. [I]) by quantum tunneling (or thermal activation) will occur 
through nucleation of hadronic bubbles in the QGP. Due to the plasma expansion, it is 
conceivable that the a field inside a bubble might get trapped at a ~ fa as the interaction 
rate of a is small compared to expansion rate. When the expansion slows down, a will start 
to oscillate and relax via production of the particles, to which the a particle is coupled, 




(1) 




(2) 



Fig. [I]. Due to the nonlinearity of U(a), the model bears nontopological soliton solutions 
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to the true vacuum. This nonequilibrium stage will take place for a time-scale of about 



therein). Typically, during this stage the dominant channel for the relaxation is the pro- 
duction of strong interacting particles such as a quanta, gluons and other mesons. These 
particles will be immediately trapped and thermalized in the nuclear matter. Perhaps the 
a quanta are short-lived to decay into photons before thermalization. If this happens, the 
decay photons would be a source for nonequilibrium photons. However, our current interest 
is the direct photon production driven by the oscillations of the a field due to parametric 
amplification. We will ignore the hydrodynamic expansion and simplify the dynamics by 
adopting a "quenched" phase transition from an initial state in local thermodynamic equi- 
librium at a temperature slightly above the critical temperature cooled instantaneously to 
zero temperature. This "quench" approximation is far from a definitive description of the 
dynamics. Nevertheless, it captures the qualitative features and allows a simple but concrete 
calculation 

At tree level the a field is inert to electromagnetic interaction, but it can couple to photon 
through a quark loop. After integrating out the quark field, we obtain an one-loop effective 
Lagrangian including strong interacting final states, 



where F^ u = d^A u — d u A^ is the photon, and G^ u = d^B v — d v B^ is a quark- ant iquark 
vector meson with mass my. It is straightforward to calculate the coupling constant g ~ 
2.6 x 10~ 3 . But the coupling constant h would depend on the wave function at the origin 
of the vector meson, \Ev(0). Assuming my ~ 2 GeV and taking a conservative value for 
|\I>y(0)| 2 ~ 0.1 fuT 3 we find h ~ 1.8 x 1CT 2 . Note that effective strong interacting 
vertices such as agg (where g denotes a gluon) and oixix should also appear at one-loop 
level. However, we have omitted them in the effective Lagrangian (|3|) because their effects 
to the particle production can be fully represented by the a strong self-couplings. Since we 
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10 fm, as revealed by hydrodynamical simulations (for example, see Ref. || and references 




(3) 



are concerned with photon production only, we integrate out the vector meson and obtain 



C = \d^a - U(a) - -F, V F^ - /-(<r - a c )F^ - Jf—d^oF^F^ (4) 
Z 4 4<7 C oTTiyCTr 

where we have dropped off higher derivative terms. At this stage, we should emphasize that 
the above obtained effective Lagrangian for the processes such as a' — *> 27 as well as 2a f — > 27 
( here a' is the shifted field with a = a'-\-a c ) is clearly understood in perturbation theory. It is 
not certain whether this effective Lagrangian can also describe the nonequilibrium situation. 
In other words, the effective vertices that account for the above mentioned processes may 
be modified in the strongly out-of-equilibrium situation. The fuller study of off-equilibrium 
effective vertices by integrating out the quark fields and vector meson is a challenging task 
that lies beyond the scope of this paper, but certainly deserves to be taken up in the near 
future. However, in this work, we will use this effective Lagrangian (|]) to compute the 
nonequilibrium photon production. 

Following the nonequilibrium closed time path formalism |5|-0, the nonequilibrium ef- 
fective Lagrangian is given by 



C = C 



c 



(5) 



where + (-) denotes the forward (backward) time branches. We then split a ± into a mean 
field and the quantum fluctuating fields: 



a ± (x,t) = (j)(t)+x ± &t), (6) 



with the tadpole conditions, 



( x ± (x,t)) = 0. (7) 



This tadpole conditions will be imposed to all orders in the corresponding expansion to 
obtain the nonequilibrium equations of motion from the first principles approach. 

To derive the nonequilibrium evolution equations that consistently take into account 
quantum fluctuation effects from the strong a self-interaction, we adopt the following Hartree 
factorization for x implemented for both ± components 



X 4 -> 6(x 2 )x 2 + constant, 
X 3 ^3( X 2 )X- 



(8) 



The expectation value will be determined self-consistently. Before proceeding any further, 
it is worth noting that although the above Hartree factorization is uncontrolled in this 
effective field theory that involves a single scalar field 0-0], our justification of using this 
approximation is based on the fact that it provides a non-perturbative framework that allows 
us to treat the strong a dynamics self-consistently. 

After doing this factorization, the Lagrangian then becomes 



C 



-C \<f>(t) + X -,A 



.(d X + ) 2 -U(t) X " 



-M 2 (t) X +2 - -F+F + ^ 



-—(Mt) - a c )F+F +lxv - -^- X + F +tlv F l 



h 2 



im v a 2 



<p 2 {t)F +0l F +0 



h 2 



<j>(t)d aX + F +0i F +a 



Am v al 



-^d, X + d aX + F + ^F" 



(9) 



where 



u{t) = 4>(t) + 



a - + ~0 2 (t) + ~£(t) 



M 2 Jt) 



W(t) + U 2 (t) + ^(t), 



m = (x 2 Kt)-(x 2 )(o). 



(10) 



Here, we have performed a subtraction of (x 2 ){t) at t = absorbing (x 2 )(0) into the finite 
renormalization of a. Besides, since we are interested in the processes of direct photon 
production driven by a time dependent (j)(t) field, in which the photons do not appear in the 
intermediate states, to avoid the gauge ambiguities, we can work on the Coulomb gauge, 
and concentrate only on physical transverse gauge fields, At [[7]. 

With the above Hartree-factorized Lagrangian in the Coulomb gauge, we perform a 
perturbative expansion in the weak couplings g and h 2 . However, the strong a dynamics 
is treated non-perturbatively JT]. Following from the tadpole conditions (|7|), we obtain the 
following full one-loop equation of motion of 4>(t) given by 
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0(0 



a- b -cf>(t) + ^ 2 (t) + ^(t) 



<f>(t) - ^(\d,A T \ 2 )(t) - 



/i 2 



0(t) + 0(t)^ 



<|I T | 2 >(t) = 0. 



(11) 



The Heisenberg field equations can be read off from the quadratic part of the Lagrangian in 
the form 



~V 2 + M 2 (t) 
l + ^-(0(t)-a c ) + 



x(t) = 0, 

h 2 



A T {t) + 



h 2 



2m 2 v a 2 



mm 



A T (t) - 
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l + il(0(O-<7 c 



V 2 l T (t) = 0. 



(12) 



Now we decompose the fields x an d A T into their Fourier mode functions U^(t) and V x %(t) 
respectively, 



X(x, t) 



2(2tt)3^ 



a n U n (t)e*"- x + h.c. 



Afe 



X %J j2{2n) 3 u Ak 



d 3 k e. 



6^(0^ + h.c. 



(13) 



where and 5 A g are destruction operators, and e A £ are linear polarization unit vectors. 
The frequencies u xk * and u Ak * can be determined from the initial states and will be specified 
below. Then the mode equations are 

d 2 



U k (t) = 0, 

h 2 



dt 2 + 



h 2 



2rriycrl 



d_ 

~dt 



+ k 2 
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1+ ^-{(j>{t) - a c ) 

Or 



-kt)+ 

(Jr. 



V Xk (t) = 0, 



(14) 



with the vacuum expectation values given by 

^ d 3 k 



(x 2 )(t) 



Mr| 2 )(t) = E 



2(27i) 3 uj xk 
A d 3 k 



U k (t)\ 2 , 



2(27i) 3 uj Ak 



\V Xk (t)\ 2 - k 2 \V Xk (t)\ : 



(15) 
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where we set the cutoff scale A ~ my, and the expectation value of the number operator for 
the asymptotic photons with momentum k is given by 

1 



<n,m> 2k 



A T (k, t) ■ A T (-k, t) +k 2 A T (k, t) ■ A T (-k, t) - 1 
= ^E[\Vx k (t)\ 2 + k 2 \V Xh (t)\*]-l. (16) 

This gives the spectral number density of the photons produced at time t, dN(t)/d 3 k. To 
solve the evolution equations flllUI3j) , we propose the following initial conditions for the 
mode functions at the time of "quench" @: 

U k {0) = 1, I7 fc (0) = -iu xk , u 2 xk = k 2 + a + bfa + ~$ ; 

V Xk (0) = 1, ^Afc(O) = -iuJAk, u A k = k, (17) 

where the initial mode functions are chosen to be at zero temperature inside the hadronic 
bubbles with the mean field <fi(t) displaced initially away from the equilibrium position (i.e., 
4>(0) = <pi 7^ 0). The above specified initial conditions are physically plausible and simple 
enough for us to investigate a quantitative description of the dynamics. 

Let us use the FL model parameters: a = 1.6 fm -2 , b = 69 fm , c = 500, and / = 9.57. 
This implies that a c = 0.36 fm -1 . This set of parameters has been used for fitting hadronic 
mass spectrum ||, and with this potential (@) (see Fig. |l|) the phase transition is weakly 
first-order at finite temperature |fL0"l . We choose the initial amplitude 0, = 0.07 fm -1 to 



solve Eq. flu]). 

In Fig. H the time evolution of the mean field <p(t) is shown. It can be seen that (f> is 
oscillating about a mean value of about 1.1 fm -1 , which is significantly different from the 
classical value a c . Also, (f)(t) oscillates with a frequency uo^ ~ 8.8 fm -1 . This is due to the 
quantum fluctuation effects coming from the a field and the gauge field At to the motion of 
4>(t) (|TTD. The back reactions of the quantum fields also account for the damping of 4>(t) with 
time. In Fig. [| we plot the time dependence of the photon number density integrated over 
momentum k, N(t). As expected, N(t) oscillates with the same frequency as <p(t). Although 
photons are produced efficiently during the first half of an oscillation, almost all of them 
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are re-absorbed to the background during the second half of the oscillation. However, it is 
important to note that the time average of N(t) is indeed increasing with time, that is to 



It is useful to calculate the time-averaged invariant photon production rate, kdR/d 3 k, 
where 



over a period from the initial time to time T. The results are shown in Fig. f| with T = 
3, 5, 10 fm. Also shown are the thermal photon production rate from a quark-gluon plasma 
and a hadron gas taken from Ref . |T^] . It is apparent that the oscillations with frequency uj^ 
produce non-thermal photons whose spectrum peaks around two photon momenta, k = 
and k = u^, which correspond to the unstable bands (or resonant bands) arising from the 
fact that the mode equations of At (|14D consistently depend on the time dependent mean 
field 4>(t) [0. Clearly, the photon production mechanism is that of parametric amplification. 
From the mode equations (0), we can easily recognize that the 4.4 fm -1 peak is resulted 
from the coupling aF 2 while the 8.8 fm -1 peak is from the interaction (da) 2 F 2 . Note that 
the 4.4 fm -1 peak has a peak value of the production rate being comparable to that of 
thermal photons, while interestingly the 8.8 fm -1 peak is almost two orders of magnitude 
larger than the thermal photons. Although these results more or less depend on the choice 
of the FL model parameters, we do not find a change of two orders of magnitude. So these 
high-energy non-thermal photons can be a distinct signature of QGP formation. 

In summary, we have computed the nonequilibrium photon production during the quark- 
hadron phase transition at high-temperature in the Friedberg-Lee model of quark-hadron 
physics. Under the "quench" approximation, the invariant production rate for nonequilib- 
rium photons driven by the oscillation of the a field due to parametric amplification is given, 
which is two orders of magnitude larger than that from a thermal quark-gluon plasma for 
photon energies around 2 GeV. These high-energy non-thermal photons can be a potential 
test of the formation of quark-gluon plasma in relativistic heavy-ion-collision experiments. 



say, photons are effectively produced from cf) oscillations due to parametric amplification || . 




(18) 
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Of course, in order to compare the results with experimental data on direct photons, a more 
realistic dynamics of the phase transition should be considered and the nonequilibrium pho- 
ton production rate should be convolved with the expansion of the plasma. This work is in 
progress. 
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of D.S.L. (K.W.N.) was supported in part by the National Science Council, ROC under the 
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FIGURE CAPTIONS 



Fig. [I]. Scalar potential U(a) with a = 1.6 fm 2 , b = 69 fm l , and c = 500, and <pi is the 
initial value. 

Fig. 0. Time evolution of the mean field <j)(t). 

Fig. Time evolution of produced total number density of nonequilibrium photons from 
quark-hadron phase transition. 

Fig. ^. Spectral production rate of nonequilibrium photons from quark-hadron phase tran- 
sition, calculated from Eq. fll8|) with T = 3,5, 10 fm, drawn with solid lines. Dashed lines 



are the rates for hadron gas and quark-gluon plasma taken from Ref. |13 . 
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FIGURES 




FIG. 1. 
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FIG. 3. 
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